Localization length of a soliton from a non-magnetic impurity in a general 

double-spin-chain model 
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A localization length of a free-spin soliton from a non-magnetic impurity is deduced in a general 
double-spin-chain model (J0-J1-J2-J3 model). We have solved a variational problem which employs 
the nearest-neighbor singlet-dimer basis. The wave function of a soliton is expressed by the Airy 
function, and the localization length (f) is found to obey a power law of the dimerization (J2 — 
J3) with an exponent —1/3; £ ~ (J2 — Js) -1 / 3 . This explains why NaV20s does not show the 
antiferromagnetic order, while CuGeC>3 does by impurity doping. When the gap exists by the bond- 
dimerization, a soliton is localized and no order is expected. Contrary, there is a possibility of the 
order when the gap is mainly due to frustration. 
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A possibility of the superconductivity upon doping car- 
riers to the spin-ladder system has attracted much in- 
terest to the quasi-one dimensional quantum spin sys- 
tem these years Jj Without the doping, the ground state 
of the spin-ladder has strong singlet-dimer correlations 
along the rung pairs. It leads to the finite energy gap and 
no long-range magnetic order. Many intensive investiga- 
tions were done in the course of doping impurities to the 
spin-gapped systems after this prediction. As its byprod- 
uct, an unexpected anomaly was found. That is, the 
long-range antiferromagnetic (AF) order appears by dop-i 
ing non-magnetic impurities to a spin-gapped material.0 
Copper sites of the base compound, SrCu2 03, is replaced 
with non-magnetic zinc atoms. The impurity destroys 
the spin gap with only 1% concentration, and causes 
the AF order at low temperature. This doping effect 
suddenly became a big topic both experimentally and 
theoreticaUyJScI The similar phenomena was also o 
served in inorganic spin-Peierls compound, CuGe03 
This compound can be explained theoretically by the 
one-dimensional frustrated spin system with a ratio of the 
next-neacest-neighbor to the nearest-neighbor interaction 
a 0.35.113 Thus, the origin of the gap is frustration. On 
the other hand, another spin-Peierls compound NaV2C>5 
has not shown the antiferromagnetic order, by the Na de- 
pletion, while the energy gap disappear.!!!! The origin of 
the gap of this system is the bond-dimerization£3 We 
can attribute this difference to the origin of the gap. 

Generally, one non-magnetic impurity atom destroys 
one singlet dimer pair, and makes one free spin- 1/2. We 
call it a soliton hereafter. Solitons should interact with 
each other to be ordered antiferromagnetically. Thus, a 
localization length of a soliton, £, should exceed a mean 
distance between impurities. Up to now, an explicit ex- 
pression for the wave function of a soliton or for its lo- 
calization length has not given analytically. Theoreti- 
cal investigations have been done mainly by numerical 
methods. cTEI In this paper, we deduce £ as a function of 
the strength of interaction bonds in a general double- 



spin-chain model by using a variation which employs the 
nearest-neighbor singlet-dimer basis. The result explains 
why a soliton localizes in the system whose energy gap 
exists by the bond-dimcrization. 

We consider the following spin Hamiltonian as the base 
system under the open boundary conditions: 
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FIG. 1. (a) Shape of the a general double spin-chain model 
with a single non-magnetic impurity at the origin, cro . (b) An 
equivalent bond-depleted model to the site-depleted model 
(a). Soliton sites are indicated by numbers, (c) A schematic 
picture for the variational basis ipi for i = — 1. An ellipse 
denotes a singlet-dimer state. 

Figure |l|(a) shows the shape of the lattice. A non- 
magnetic impurity is located at the <To-site. This site- 
depleted lattice is equivalent to the bond-depleted lat- 
tice as shown in Fig. 0(b), if we shift spins on the left 
side of the impurity by one lattice spacing to the right 
and exchange the upper and the lower chain. We make 
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a variational analysis to this bond-depleted model. 

The ground state of the base system has a strong 
nearest-neighbor singlet-dimer correlation when it has 
a finite energy gap. We consider this picture remains 
valid even after an impurity is doped, i.e., the ground 
state can be described in terms of one soliton and 
27V nearest-neighbor singlet-dimer states. Far-neighbor 
dimer states are neglected. In this sense, the analysis 
is variational. This soliton approach has already suc- 
ceeded in explaining propetties-pf the low-lying states in 
the fully-frustrated system. EjEj The present approxima- 
tion is good as long as the nearest-neighbor interactions, 
J 2 and J3 dominate the other bonds. 

The variational basis is as depicted in Fig. [j] (c) and 
is written as follows. 

Ipi = [c-N, T-n] ■ ■ ■ [n-l,(Ti] Tt fo+l, <7i+l] ■ ■ ■ b~N, Cjv], 

where i > 0, and [a, b] = (Talfc — laT&)/v2 denotes the 
singlet-dimer state. We convention that the spin-a is al- 
ways located in the upper chain, and the spin-& in the 
lower chain. The basis ipi for i < is defined in the 
same way. A location of a soliton | j is restricted within 
the upper chain, thus there are (2N + 1) different states. 
Note that a state with a soliton in the lower chain can 
be expanded in a linear combination of states with a soli- 
ton in the upper chain. The present variational basis is 
not orthogonal to each other and satisfies the following 
relation: 
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The matrix elements of the Hamiltonian, (i/)i\'H\ipj) 
different forms in regard to a sign of i and j : 
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take 



= -max(i, j)(J 2 - J3) + j\{J 2 - J\- Jo) 
+2NJ 2 + J Q (l-S iJ ), 

for max(i,j) > and min(i, j) > 0, 
= \i- j\{Ja -Ji- Jo) + 2NJ 2 -J 3 + 2(Jx + 

for max(i,j) > and min(i,j) < 0, 
= -min(«, j)(J 3 -Ji- Jo) + 2NJ 2 + Jo + 2Ji 

for max(i,j) = and min(i,j) < 0, 
= max(i, j)(J 2 - J 3 ) + \i-j\(J3 - Ji - Jo) 
+2NJ 2 + 2J 1 -J 3 -J 1 S itj , 

for max(«,j) < and min(i, j) < 0, 
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where max(i, j) = i for i > j, min(i,j) = i for i < j, and 
5jj is the Kronecker delta. The constant energy term is 
lower in the case of i,j > than that in i,j < 0, when 
Ji < J3. This makes a soliton confined within the region 
of positive sites in the low-energy limit. Therefore, an 
impurity acts as a potential barrier to a soliton. 

Our task is to find a function, \t var = X^jCiV'ij that 
minimizes the energy expectation, 
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If we diagonalize the denominator and rewrite the nu- 
merator with its eigenfunction, this variational problem 
is transformed into a simple eigenvalue problem. We first 
solve this transformed eigenvalue problem by the numer- 
ical diagonalization for a finite system with N = 200. 
After that, we give an analytic solution in the continu- 
ous limit, N — > 00, when the wave number of the lowest 
energy state is zero. 
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FIG. 2. The ground-state wave function obtained by 
numerical diagonalization of the variational problem for 
(a): Jo = 0.35, (b):J = 0.5, and (c):J = 0.7. J 3 is 
set unity and Ji = Jo. The strength of the dimerization 
J2 - J3 = 2.441 x 10" 4 . Solid lines in (a) and (b) are the 
asymptotic form of the Airy function. 
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Figure || shows behaviors of the ground-state wave 
function. We set J3 unity and vary Jo = Ji as (a): 
Jo = 0.35, which corresponds to CuGeC>3; (b): J = 0.5, 
the Majumdar-Ghosh model; and (c): Jo = 0.7, which is 
in the incommensurate phase. Since our analysis is base 
on the nearest-neighbor singlet-dimer state, the results 
in the incommensurate phase is not beyond the spec- 
ulation level. For Jo < 0.5, the wave function of the 
non-dimerized case (J2 = J3) is the sine- function with 
a phase shift caused by the impurity. The phase shift 
decreases with a decrease of J\, since the effective po- 
tential barrier at the impurity increases, and it becomes 
an open edge at J\ — 0. A soliton can be regarded as 
a free particle in a potential well in this case. When 
the dimerization is switched on, the wave function sud- 
denly collapse to the impurity and exhibits an exponen- 
tial decay. The strength of the dimerization, J2 — J3, 
is only 2.441 x 10~ 4 in this figure. Therefore, a soliton 
becomes localized even with a sufficiently weak dimeriza- 
tion. It can be noticed that the wave function is more 
localized in the case Jo = 0.5 as compared with that of 
Jo = 0.35. We will clarify this Jo-dependence of £ in 
the analytic solution. In the case of Jo = 0.7, the wave 
function exhibits an oscillation with an incommensurate 
wave number, k[ nc , and is already localized without the 
dimerization. The dimerization effect only decreases a 
fraction of the localization length. Since the origin of the 
localization in Jo < 0.5 can be attributed to the bond- 
dimerization alone, the localization phenomenon in the 
incommensurate phase is considered quite different from 
the former ones. We leave to clarify this problem for the 
future investigation. 

The variational problem can be solved analytically by 
the Fourier transformation, since the denominator is di- 
agonalizcd in the N — > 00 limit. Therefore, we rewrite 
the problem with a new basis, \4>k) = X]„ ex P[^ n ]l^n)- 
The basis relation becomes 



5 — 4 cos k 



S k> i + 0(1/N). 
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The leading and the second terms of the diagonal part of 
the Hamiltonian, TLk,k = ( < Pk\'H\^k), are 
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5 — 4 cos k 



+ (J + Ji) + 8(J 2 - Ji - J ) 



4 — 5 cos k 
(5 - 4 cos A;) 2 
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where we have dropped the constant term, 2NJ 2 - The 
leading term of the off-diagonal part, TLk,ii is 



3 (J 2 - J 3 ) 

4 2N + 1 

3 ( 1 
2 



AN 



5 cos((fc + l)N) - 4 cos((fc - l)N) 
(5-4 cos k) (5 - 4 cos I) 
1 \ 1 - cos((fc - l)N)~ 



4 cos k 5 — 4 cos I 



1 — cos(fc — I) 



The off-diagonal parts are always smaller than the diag- 
onal parts by an order of N. 



Before proceeding to the solution, we list several no- 
tices. First, the wave number of the ground state must be 
zero, because we consider the continuous limit, N — > 00 
and k — > 0. In the case of Jo = Ji, and J 2 = J3, the 
wave number that minimizes the variational energy, fci nc , 
is zero for Ji/J 2 < 9/17, and otherwise it is given by 



cos kh 



\ (5 - V9(2 Ji - J2)M) . (11) 



This incommensurate point was^aJso obtained by the vari- 
ational matrix-product ansatz$3 Secondly, we find the 
diagonal part (||) diverge for J 2 — J3 7^ 0. In order to 
avoid this divergence, we introduce a cutoff factor to the 
momentum as k — > fe + i8, which has no effect to the 
physical results by taking a limit S — > after N — > 00. 

We rewrite the matrix elements and pick up only the 
leading term of the off-diagonal part and terms up to k 2 
in the diagonal part. Then, we have the continuous limit 
of the Hamiltonian TLk 1 ■ 



7~Lk,i = 



const. + 

3 J 2 - J 3 

4 2N+ 1 
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1 
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where the constant term is (3J 3 + 6J + 3Ji — 8J 2 )/8. 
Apart from this constant term, the Hamiltonian is equiv- 
alent to the following in the real space representation: 



2m dx 2 



+ 7W2- J 3 )|z|exp[-<5|x|]. (13) 



Here, x is the distance from the impurity, and m is an 
effective mass of a soliton given by 



- 1 = (J + Ji)/2 + 9(J 2 - J x - Jo). 



(14) 



We take the limit 8 — > at this stage. The first term 



of the eq. (13) is the kinetic energy of a soliton, and 
the second term is the triangular potential which comes 
from the bond-dimerization. By a scale transformation, 
X = 9x, the eigenvalue equation = Eq^ becomes 



dX 2 



(X-E') 



* = 0, 



(15) 



where E' = E c x 2md~ 2 with 6» 3 = 3m(J 2 - J 3 )/2. A 
solution of this equation is known as the Airy function- 
* = Ai(X - E'), with the first eigenvalue E' = 2.33813 
Then, the energy, Ec, and the localization length, £, of 
a soliton is estimated as 



Ec = 



E' 



1.532 x m- 1/3 (J 2 - J 3 ) 2/3 , (16) 
( 10 ) £~(3 + .E')x0~ 1 = 4 -663xm- 1 / 3 (J2- ^) _1/3 - (17) 
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The numeric factor, (3+E'), is determined since the Airy 
function decays exponentially for large X as 



3 



Mx) ~ ^v^ exp[ " 2x3/2/3] ' (18) 

and Ai(3) ~ 0.01. We plot this asymptotic form in Fig. 
^|(a) and (b) and compare with the numerical solution. 
The agreements are excellent, even though we have ne- 
glected the impurity term of the off-diagonal matrix ele- 
ments in the analytic solution. As depicted in Fig. |l| (b), 
the Jo- and the Ji-bond vanish and the J2-bond changes 
to the Ji at the impurity site. This means the effective 
mass given by eq. ( |l4|) is reduced to mT 1 = 9J\ only at 
this site. Therefore, the impurity attracts a soliton. The 
analytic localization length, eq. (0), should be modified 
as £ — a by a some constant term a regarding the phase 
shift. We will determine a by fitting £ — a to the nu- 
merical results of the variation in the following, since the 
phase shift can not be obtained by the present analysis. 

Figure || compares the analytic results of £ to those 
obtained by the numerical diagonalization. The numer- 
ical data are estimated at the point where the absolute 
value of the normalized wave function l^varCOl becomes 
smaller than 0.01. Analytic data depicted by lines are 
£ — a with a = 5.5. This small subtraction merely im- 
proves the fit in a small-£ region, £ < 20. It does not in- 
fluence the asymptotic form of £ which is in a power law 
with an exponent —1/3. We have checked this behavior 
remains in a wide region of < J0/J3 < 0.5. We can also 
predict the strength of the dimerization, S = (J2 — </3)/2, 
in CuGeC>3 by a plot of Jo = 0.35 of this figure. Since 
the AF order appears at the 2% doping but disappears 
at the 1%J3 £ would be 12.5 < £ < 25, which roughly 
corresponds to 0.003 < 5 < 0.025. 
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FIG. 3. Localization length estimated from the wave func- 
tion obtained by numerical diagonalization for Jo — J\ — 0.35 
(o) and Jo = Ji = 0.5 (x). J3 is set unity. Lines are the an- 
alytic results of £ — a with a = 5.5. 

In summary, we have given an explicit expression for 
the localization length of a soliton when a non-magnetic 
impurity is doped in a general double-spin-chain model. 
The expression (|l7|) classifies the behavior of a soliton 
by the origin of the energy gap. When the gap origi- 
nates in the bond-dimerization, a soliton is strongly lo- 
calized near the impurity, which means there expects no 



antiferromagnetic order. On the other hand, a soliton 
can interact with another soliton when the gap is formed 
mainly by frustration, since £ diverges. This result ex- 
plains why NaV2C>5 does not exhibit the AF order by 
the Na depletion, while CuGc03 does by the Zn doping. 
The present analysis clarified that the motion of a soli- 
ton is well-described by a quantum particle in a potential 
well. A soliton is confined to one side of the impurity by 
the potential barrier proportional to J3 — J\, and thus, 
it can penetrate to the other side as J\ approaches J3. It 
also feels a triangular potential which is an outcome of 
the bond-dimerization. We conclude that this triangular 
potential is the cause of the localization of a soliton in 
the dimer phase. Another type of the localization was 
also found in the incommensurate phase. It occurs ir- 
respective of the bond-dimerization. Derivation of the 
interaction between solitons and the localization prob- 
lem in the incommensurate phase are left for the future 
investigations. 
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